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Abstract: The behavior of the structure function at high energies 
(high densities) is directly associated to the gluon distribution. In this paper 
we analyze the asymptotic behavior of the structure function at very high 
densities considering a leading and a higher twist relation between these 
two quantities. We verify that (a) if the leading twist relation between the 
structure function and the gluon distribution is used, we recover the black 
disc limit, and (b) a softer behavior is predicted if a higher twist relation 
between F2 and xG is considered. While the first behavior is well-established 
by the current high density approaches, the higher twist behavior is a new 
result. In both cases, the F2 structure function unitarizes and the Froissart 
boundary is not violated in the asymptotic regime of high density QCD. 
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1 Introduction 



The physics of high-density QCD has become an increasingly active sub- 
ject of research, both from experimental and theoretical points of view. In 
deep inelastic scattering (DIS) the parton high density regime corresponds 
to the small x region and represents the challenge of studying the interface 
between the perturbative and nonperturbative QCD, with the peculiar fea- 
ture that this transition is taken in a kinematical region where the strong 
coupling constant ag is small. By the domain of perturbative QCD we mean 
the region where the parton picture has been developed and the separation 
between the short and long distance contributions (the coUinear factoriza- 
tion) is made possible by the use of the operator product expansion (OPE). 
The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations are 
the evolution equations in this kinematical region. These equations are valid 
at leading twist, i.e. at a large value of the photon virtuality Q^, where a 
subclass of all possible Feynmann graphs are dominant and the hard coeffi- 
cient function is connected to the proton by only two parton lines (For more 
details see e.g. Ref. 0]). For small values of Q^, this picture has corrections 
predicted by the OPE that contribute at relative order O^l/Q"^) and beyond 
[0{-^)"', n = 2,3, ...]. These are commonly called higher twist corrections. 

In the limit of small values of x (< 10~^), on the other hand, one expects 
to see new features inside the nucleon: the density of gluons and quarks 
become very high and an associated new dynamical effect is expected to stop 
the further growth of the structure functions. In particular, for a fixed hard 
scale ^ ^qcd^ ^^e OPE eventually breaks down at sufficiently small x [0]. 
Ultimately, the physics in the region of high parton densities will be described 
by nonperturbative methods, which is still waiting for a satisfactory solution 
in QCD. However, the transition from the moderate x region towards the 
small X limit may possibly be accessible in perturbation theory, and, hence, 
allows us to test the ideas about the onset of nonperturbative dynamics. 

At this moment, there are several approaches in the literature that pro- 
pose distinct evolution equations for the description of the gluon distribution 
in high density hmit 0, |^ B H • general these evolution equations resum a 
class of higher twist diagrams dominant in this kinematical region that con- 
tain powers of the function k(x,Q'^) = ^'^^^i^ ^^^^^f \ which represents the 
probability of gluon-gluon interaction inside the parton cascade. Moreover, 
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these equations match (a) the DLA hmit of the DGLAP evolution equation 
in the hmit of low parton densities (k — > 0); (b)the GLR equation and the 
Glauber-Mueller formula as first terms of the high density effects. The main 
differences between these approaches occurs in the limit of very large den- 
sities, where all powers of k should be resummed. Although the complete 
demonstration of the equivalence between these formulations in the region 
of large k is still an open question, some steps in this direction were given 
recently p, |10 . 



Our goal in this paper is not to demonstrate the equivalence between the 
distinct approaches, but to analyze the relation between the gluon distribu- 
tion xG and the F2 structure function in the asymptotic regime of very high 
densities (energies). In general, a leading twist relation is used to relate these 
quantities and the solutions of the evolution equations used as input in the 
calculations. However, as discussed above, the high density approaches re- 
sum a class of higher twist diagrams, which implies that we should careful in 
our predictions for the quantities which would be measured. Here we study 
the behavior of F2, predicted by a leading and a higher twist relation, at 
large densities, where analytical calculations are possible. In this regime we 
expect the blackness of the cross section, as predicted by Gribov many years 
ago |jll|]. We verify that this behavior is recovered by the current high density 



approaches, if a leading twist relation between F2 and xG is used. However, 
a softer behavior is obtained when we consider that this relation is modified 
by the higher twist terms associated to the high density, demonstrating the 
importance of this contribution. In both cases, the Froissart boundary is not 
violated. 

The paper is organized as follows: in the next section we briefly review 
the derivation of the leading and higher twist relations between F2 and xG, 
as well as the black disc limit. After, in section 3, we derive the asymptotic 
solution of the AGL equation at fixed and running and compare this 
solution with the predictions of the McLerran - Venugopalan high density 
approach. Finally, in section 4, we present our results and conclusions. 
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2 The F2 Structure Function and the Black 
Disc Limit 



We start from the space-time picture of the ep processes ||rT|. The deep 
inelastic scattering ep —>■ e+X is characterized by a large electron energy loss 
u (in the target rest frame) and an invariant momentum transfer = — 
between the incoming and outgoing electron such that x = Q'^ /2mMV is fixed. 
In terms of Fock states we then view the ep scattering as follows: the electron 
emits a photon (|e >— |e7 >) with = u and p'f^ ~ Q^, after the photon 
splits into a qq (|e7 >— \eqq >) and typically travels a distance Ic ~ 1/mj^x, 
referred as the coherence length, before interacting in the nucleon. For small 
X (large s, where ^/s is 7*p center-of-mass energy), the photon converts to a 
quark pair at a large distance before it interacts to the target. Consequently, 
the space-time picture of the DIS in the target rest frame can be viewed as the 
decay of the virtual photon at high energy (small x) into a quark-antiquark 
pair long before the interaction with the target. The qq pair subsequently 
interacts with the target. In the small x region, where x <C 2mR^ 1^ 
pair crosses the target with fixed transverse distance between the quarks. 



Following Gribov |]TT|, we may write a double dispersion relation for the 
forward 7*p elastic amplitude A, related to the total cross section by the 
optical theorem {ImA = sa{s, Q^)), and obtain for fixed s 

-is,Q')=j:J Q. Pis, M\ M-)\lm A,,,,is, M^ M-) , (] 

where M and M' are the invariant masses of the incoming and outgoing qq 
pair. If we assume that the forward qq + p scattering does not change the 
momentum of the quarks then Agg+p is proportional to S{M^ — M'^), and (|^) 
becomes 

where the spectral function p(s, M^) is the density of qq states, which may be 
expressed in terms of the 7* — >■ qq matrix element |jl2||. Using that = {k^ + 



m^^/[z{l — zy\, where kt and z are the transverse and longitudinal momentum 
components of the quark with mass mg, we can express the integral over the 
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mass M of the qq in terms of a two-dimensional integral over z and kt- Instead 
of kti it is useful to work with the transverse coordinate (impact parameter 
representation), which is the variable Fourier conjugate to kt, resulting [IT3 



F2(x,Q2) = -^a{s,Q') = fdz f -^\^(z,r,)\' a,,{z,n) , (3) 

where 

l^(^,n)r = ^i:e?{[-^ + (1 - zf\^K,{.nf^m\K,{enf} . (4) 

The photon wave function \E'(z,rf) is simply the Fourier transform of the 
matrix element for the transition 7* — > qq. Moreover, is the electro- 
magnetic coupling constant, = z{\ — z)Q'^ + mf, rrii is the quark mass, 
Uf is the number of active flavors, is the square of the parton charge (in 
units of e), Kq i are the modified Bessel functions and z is the fraction of 
the photon's light-cone momentum carried by one of the quarks of the pair. 
In the leading log(l/x) approximation we can neglect the change of z during 
the interaction and describe the cross section (j''^(2;, 4/r^) as a function of the 
variable x. Considering only light quarks {i = u, d, s) F2 can be expressed 



by [14 



1 V- 2 /"q? (^^^t 

u,d,s 

Using that ag-g = g(3a,(^)/4) vr^ ^^(a;, ^) IH, where xG{x,^) is the 



Ca _ 
nucleon gluon distribution, we get 

The above expression explicit the direct dependence of the F2 structure func- 
tion with the behavior of the gluon distribution, which is related to the QCD 
dynamics at high energies. In the linear regime, where the high densities can 
be disregarded, the small x behavior of the gluon distribution is given by 
the solutions of the DGLAP [jl| and/or BFKL evolution equations. The 
common feature of these equations is the steep increase of xG as x decreases. 



4 



This steep increase cannot persist down to arbitrary low values of x since it 
violates a fundamental principle of quantum theory, i. e. the unitarity. In the 
context of relativistic quantum field theory of the strong interactions, uni- 
tarity implies that the cross section cannot increase with increasing energy 
s above log^ s: the Froissart's theorem |T^. In the next section we will con- 



sider the modifications in the behavior of the gluon distribution associated 
to the high density effects. The presence of a high parton density implies a 
slow growth of xG at high energies, and consequently from Eq. (^), of the 
structure function. However, the Eq. (P) is a leading twist relation which 
will eventually breaks down when we consider the higher twist terms in the 
evolution. 

We now consider that the relation between the structure function and the 
gluon distribution is modified by the particular type of higher twist terms 
associated with the high density corrections (See also [T^ for a similar calcu- 
lation in the infinite momentum frame). In general, the higher twist contri- 
butions should be significant at small x and Q^, implying large perturbative 
corrections to the conventional leading twist relations in this kinematical 
region. We estimated these corrections considering the s-channel unitarity 
constraint in the interaction cross section of the quark-antiquark pair with 
In this case the structure function is given by 



the target 117 



2^ ^ ' 

f=u,d,s 



f 



(frt 



d\{l 



aii{x,4,/rf)S{bt) 



} 



(7) 



A similar expression was used in Ref. [|18[ for a phenomenological analysis of 
the ep process, where the geometrical structure of the collision (the bt depen- 
dence) was disregarded and xg oc a;^^ (A > 0) was assumed, resulting a very 
good description of the HERA data. The main point of the above expression 
is that this resums a large class of higher twist terms, as demonstrated in 
Ref. 
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The use of a Gaussian parameterization for the nucleon profile function 

- ''^ 

S{bt) = ^! where R is the spatial gluon distribution inside the proton 

We obtain that the Fo structure function 



{C + ln{K,{x,rl)) + E,{K,{x,r1))} , (8) 



simplifies the calculations, 
can be written as 



2^ ^ 



u,d,s 



4 
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where Hq = i/dhzc = {2as/3R'^) n rf xGisfix, jj). This equation allows to 
estimate the high density corrections to the structure function in the DLA 
limit. Expanding the equation for small Kg, the first term (Born term) 
will correspond to the usual DGLAP equation in the small x region [Eq. (^]. 

Before we discuss the dynamics for the gluon distribution at high densi- 
ties, it is important to consider a general property of the structure function 
(i.e. the cross section): the black disc limit. As shown before, when the 
coherence length Ic considerably exceeds the diameter of a target which is 
at rest, the virtual photon transforms into hadron components well before 
the target, which implies that the small x physics probes the interaction of 
various hadron wave packets with a target. The geometrical limit for the 
cross section of virtual photon scattering off a nucleon target follows from 
the assumption that a target is black for the dominant hadron components 
in the wave function of the virtual photon. An identical assumption is of- 
ten used to deduct the Froissart limit of hadron-hadron interactions at high 
energies. 

Under the assumption that the interaction is black, we can estimate the 
total cross section (H) using the following approximations: p(s, M^) oc 
and (jgg+p(s, M^) = ttR"^. Consequently, we can get a{s, M^) oc R"^ ^n[{M^^^ + 
Q'^)l [Ml^-^ + Q^)]. Using that M^.„ 4m^. whereas M^^^^ oc s, results that 

F2(x,g2) ocg'i^'/n^ , (9) 

which is the black disc limit for the F2 structure function. This bound, 
derived from a geometrical analysis, represents the maximum value allowed 
for the cross sections at high energies (high densities). In the next sections we 
will analyze the predictions of the high density approaches for the asymptotic 
behavior of F2 and verify if this bound is saturated or a softer behavior is 
predicted. 

3 The AGL High Density Approach 

About seventeen years ago, Gribov, Levin, and Ryskin (GLR) |^ performed 
a detailed study of the high density limit of QCD in the double logarithmic 
approximation (DLA). They argued that the physical processes of interaction 
and recombination of partons become important in the parton cascade at a 
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large value of the parton density, and that these high density corrections 
could be expressed in a new evolution equation - the GLR equation. Some 
years ago an eikonal approach to the high density corrections was proposed in 
the literature by Ayala, Gay Ducati, and Levin (AGL) |1^ which improves 
the GLR approach. This approach address the high density regime from the 
region where perturbative QCD is valid by summing corrections to the linear 
DGLAP evolution equation. 

As we are only interested in the solutions of the AGL equation, we refer 
the original papers [l^ for details in its derivation. Here we only present 
the AGL equation, which is given by 

^'""fyQ^^ = ^^{C + ln[KGix, Q')] + E^[KGix, Q')]} , (10) 

where C is the Euler constant, Ei is the exponential function and the function 
Kg is defined by 

^g{x,Q')^^^^xG{x,Q') , (11) 

and represents the probability of gluon-gluon interaction inside the parton 
cascade. Using the above definition for kg we can rewrite the expression (|T0|) 
in a more convenient form (for fixed as) 

d^kG{y,e) dkG{y,e) asN^ . . , 2m , 77. r ( n■^\^\ 

oyot ay n 

^ F{kg). (12) 



Analyzing the structure of the Eq. (|12]) we see that it has a solution which 
depends only on y. In it was shown that this solution is the asymptotic 
solution of the AGL equation. In this case we have that at large values of 
densities the asymptotic solution is given by 

fcc^d/) = ^y. (13) 

TT 



Substituting the definition of kg [Eq. ([TT|) ] in the above solution, in the 
asymptotic regime the behavior of the gluon distribution is given by (at 
fixed as) 

xGix,Q') = ^-^^^lnik. (14) 



7 



Therefore, the gluon distribution does not saturate at small values of x, but 
is linearly proportional to Ins (s ~ 1/x). However, this behavior is softer 



than predicted by the DGLAP equation {xG oc exp[^J Inl / x]) and the BFKL 
equation {xG oc x~^, A > 0). We obtain that the gluon distribution presents 
a partial saturation in its behavior. 

The analyzes above was made at leading order, considering a fixed cou- 
pling constant as- One question is the possible modifications in the asymp- 
totic behavior associated to next-to-leading order corrections as, for instance, 
the running of a^. In order to answer this question we need to extend the 
Refs. [P, 17 1 for running coupling constant, as we will do in the following. 



Since the QCD coupling constant is given by as{Q'^) = 47r/(/?oe), where 
Po = 11 — 2/3n/ {uf is the number of flavors) and e = /nQ^/Ag^^, the 
relation between the gluon distribution and the function kg can be expressed 
by 

xGix,Q') = ^^^eKG{x,Q') . (15) 
iVcVr 47r 

Using this result in Eq. (^) we obtain that the AGL equation with running 
ttg is given by 

+ = —;^{C + Hna{x,Q)] 

+ E,[kg{x,Q^)]} 

= H{kg) . (16) 

The AGL equation at fixed as is a direct consequence of the above equation 
in the limit of large e. Moreover, this equation also has a solution which 
depends only on y. In this case we have 

with the solution 



(y-yo). (18) 
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Consequently, the asymptotic behavior of the gluon distribution is given by 



xG{x, Q ) = ^ (-) • (19) 

As expected at large values of e (Q^) this solution reduces to the solution at 
fixed Us- The main difference is the behavior at small values of e, where the 
prefactor in (|19]) is important. However, the partial saturation of the gluon 
distribution is not modified by the running of as- It agrees with the result 
that the unitarity corrections are expected to be relevant before the next to 
leading order corrections [If. 

Before we present our results and conclusions in the next section, let us 
discuss the generality of the above results for the gluon distribution by the 
comparison with the predictions of other high density approach. In contrast 
to the AGL formalism, the approach proposed by McLerran, Venugopalan 
and collaborators (MV-JKLW approach) [^, |^ address the high density 
regime from the nonperturbative region by developing an effective Lagrangian 
approach to a very dense system. In Ref. ||2^ the authors have considered the 



DLA limit of the functional evolution equation derived from this approach, 
obtaining an equation similar, but not identical, to the AGL equation (we 
refer for details). One of the main results of this work was the solution 
of the proposed equation for large values of densities (large n) , which is given 

by 

:,G{x, Q') = ^'^^^^ ~ ^^"^ ln(-] . (20) 

2 \xj 

The remarkable feature of this solution is the same and x dependence of 
the AGL equation in the asymptotic regime. The difference in the prefactors 
is a function of the distinct normalizations and approximations used in the 
two approaches. 



4 The Asymptotic Behavior of F2 

The main conclusion of our analysis in the previous section was the uni- 
versality of the X and dependence of the gluon distribution (the partial 
saturation) in the asymptotic regime of the high density QCD. In this section 
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we will consider the consequences of the partial saturation in the behavior of 
the structure function, which is experimentally measured. 

Using the solution of the AGL equation in the asymptotic regime [Eq. 
(iMl)] as input in the Eq. (Bl) we get 



F2{x,Q')^^R'Q'ln(-) . 



(21) 



We see that the partial saturation of the gluon distribution implies that the 
F2 structure function does not saturate at small values of x, but is linearly 
proportional to Ins. Basically, we verify that the AGL approach, as well 
as the MV-JKLW approach maximizes the prediction for the F2 structure 
function in the asymptotic regime of very high densities, resulting in the 
black disc prediction, when the leading twist relation [Eq. (|^)] is used. 

We now consider that the relation between the structure function and the 
gluon distribution is given by the Eq. (H). In the asymptotic regime (large 



Kq) we obtain 



27r2 



u,d,s 



(frt 



7rr7 



ln{Kq{x,r1)) 



(22) 



Using the asymptotic solution of the AGL equation we can determine Kg 
4/9kg' at large values of densities, and so 



F2{x,Q' 



27r2 



u,d,s 



" ^In 
1 nrf 



In 



,1, 



X 



(23) 



Therefore, considering the contribution of the higher twist terms in the re- 
lation between the structure function and the gluon distribution, we predict 
the following asymptotic behavior for the structure function 



F2{x,Q' 



37r2 



In 



Aas . , 1 . 
In — 

3 ^x' 



(24) 



We see that the inclusion of the higher twist terms implies a softer dependence 
of F2 with the energy than obtained using the leading twist relation. In con- 
trast with the leading twist prediction, in this case the black disc limit is not 
maximized. However, in both cases the structure function does not violate 
the Froissart boundary in the asymptotic regime of high density QCD. The 
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demonstration of the behavior ( p^ using the other approaches for hdQCD 
is an important open question. In particular, the analyzes of the predictions 
for F2 obtained in Ref. |TB[ in the infinite momentum frame for the very 



high density regime would be a important check of our results. Our study 
demonstrates the importance of the correct connection (leading or higher 
twist relation) between the gluon distribution and the observables when we 
are analyzing the interface between perturbative and nonperturbative QCD. 

From the results for the structure function in the asymptotic regime we 
can see that this regime is characterized by the identity 

'-^^-F.M,Q-'), (25) 

which is an important signature of the asymptotic regime of high density 
QCD. This regime should be reached for the case of an interaction with 
nuclei at smaller parton densities than in a nucleon, since = As x kn, 
where ktv is given by the expression (|ll|). 



5 Conclusions 

In this paper we have analyzed the description of the F2 structure function 
in the regime of very high densities, assuming first that the leading twist 
relation between F2 and xG is valid and second that this relation is modified 
by the higher twist terms associated to the high density corrections. In the 
first case we have obtained that the corresponding F2 structure function is 
linearly proportional to Q^B? Inl/ x, which agrees with the results predicted 
in the black disc limit. This result is not high density approach dependent. In 
the second case a softer behavior is obtained. In both cases, the F2 structure 
function unitarizes and, as expected by construction, the black disc boundary 
is not violated in the asymptotic regime of high density QCD. We verify that 
when the parton density is such that the proton becomes black and the 
interaction probability is unity, the structure function becomes proportional 
to and a soft energy dependence is predicted. As a by product, we have 
derived the AGL equation for running as and obtained a signature for the 
asymptotic regime of QCD. 

Our main conclusion is that we should be careful to relate the gluon dis- 
tribution with the observables (for instance, F2, ^2^^, heavy quark production. 
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...) in the high density hmit, since besides the behavior of the gluon distri- 
bution, also the relation between these quantities must be modified from the 
usual leading twist expressions. 
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